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1. Introduction 

The old problem of magnetic structure of rare-earth metals and their compounds 
is still a subject of experimental and theoretical investigations. These substances 
have complicated phase diagrams and demonstrate a number of orientational phase 
transitions. In particular, such transitions take place in the orthoferrites and practically 
important intermetallic compounds RC05 (R=Pr,Nd,Tb,Dy,Ho), see, e.g., Ref. 
Qualitative explanation of these transitions has been obtained many years ago within 
the Heisenberg model with inclusion of magnetic anisotropy 0. In a number of 
systems, lattice (magnetoelastic) effects are important. The standard description is 
usually performed within mean-field approaches. However, quantitative comparison 
with experimental data requires a more detailed treatment. 

Provided that the orientational transition temperature is low (in comparison with 
the magnetic ordering point), spin- wave theory is applicable 0. In the simplest case of 
the second-order anisotropy the magnetization lies either along the easy axis or in the 
easy plane. Inclusion of higher-order anisotropy constants can lead to cone phases where 
magnetization makes the angle 9 with the z-axis. The case 9 = 71/ 2 was considered in 
Refs. ||, |], H || where the temperature renormalization of anisotropy constants and the 
spin-wave spectrum in Tb and Dy within the standard spin- wave theory were calculated. 

In the present paper we consider the cone phase with arbitrary < 9 < tt/2. The 
situation here is analogous to that for the field-induced orientational phase transitions, 
e.g., in the transverse- field Ising model (see Ref.[|7| and references therein). Unlike 
the latter model, one can expect that at low enough temperatures and small value 
of anisotropy the spin-wave theory is applicable for an arbitrary relation between 
anisotropy parameters, not only close to the orientational phase transition. Even for 
the second-order easy-plane anisotropy, the Holstein-Primakoff representation for spin 
operators used in Refs. || |5], |] leads to so-called kinematical inconsistencies because 
of incorrect treating on-site kinematical relations. To avoid this difficulty, we use the 
technique of spin coherent states. Our approach is to some extent similar to the operator 
approach used in Ref. J4|, but gives a possibility to treat more simply higher-order 
anisotropy constants, as well as to calculate higher-order terms of 1/ ^-expansion. 

The anisotropic Heisenberg model used is formulated in Sect 2. In Sect 3 we develop 
a special form of the l/S'-expansion which gives a possibility to take into account exactly 
on-site kinematical relations. In Sect 4 we treat the cone-plane transition owing to the 
temperature dependence of the cone angle 9 and discuss experimental data on the rare 
earth metals. 
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2. The model and mean-field approximation 

We start from the Hamiltonian of the Heisenberg model with inclusion of single-site 
anisotropy 

w = -i E + Bi E(o 2 °), + Bi Y,{o% (1) 

where J > is the exchange parameter, 

0° = 3(5 Z ) 2 - 5(5 + 1) 

0° = 35(5 2 ) 4 - 305(5 + l){S z f + 25(5 2 ) 2 

+ 35 2 (5 + l) 2 -65(5 + 1) (2) 

are the irreducible tensor operators of second and fourth orders, B™ are the 
corresponding anisotropy constants. 

Up to unimportant constant we can rewrite the Hamiltonian (|l|) in the form 

w = -i E + D E(^ 2 ) 2 + ^ E(^) 4 (3) 

(tf) i i 

where 

D = W° 2 - [305(5 + 1) - 25}B° 4 (4) 
D' = 355° 

For D,D' > spins of magnetic ions lie in the easy plane xy, while for D, D' < we 
have the easy axis z. For D > 0, D' < a first-order transition takes place between 
the easy plane (which is favored by second-order anisotropy) and easy axis (which is 
favoured by large \D'\). We consider only the case D < 0, D' > where the cone phase 
occurs at intermediate values of D/(2D'S 2 ), so that spin orientation direction makes 
the angle 9 with the z-axis and the orientational phase transitions are of the second 
order. This is the case for Gd and also for Ho, Er in low-temperature phases. 
In the phenomenological approach it is supposed (see, e.g., Refs. [p], ||) 

F = F is + D{T){S cos 6) 2 + D'{T){S cos6) 4 (5) 

where F^ s is the isotropic (^-independent) part of the free energy. Then we obtain by 
minimization of F 

so that at the point where D(T) = the spins become directed in the xy plane while 
at \D(T)\ > 2D'{T)S 2 they are aligned along the z-axis. The temperature dependence 
of D{T) is supposed to have the form 

D(T) = 2D'S 2 (T 1 - T)/(T 2 - T x ) (7) 
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with D'(T) > 0. Thus at T — T\ the transition from the easy-plane to cone structure 
takes place, while at T = T 2 the transition from the cone to easy-axis structure occurs. 
At the same time, Zener's || resut for the temperature dependence of anisotropy 
constants in an axially symmetric state with 9 = has the form 

B?(T) = B?M l W 2 (8) 

where M = (S z )/S is the relative magnetization, and D(T), D'(T) are determined 
by the same relations (§ with Bf -> Bf(T). As pointed in Refs. §, |, §, §, the 
temperature dependences of anisotropy constants have a more complicated form for 
the cone structures with 9 > (in fact, only the case 9 = n/2 is discussed in Refs. 
I SSD- 

A systematic way of calculating temperature dependences of anisotropy constants 
is the 1 / S-expansion which is considered in the next section. 



3. The 1/S expansion of the partition function 



The l/S'-expansion developed here is slightly different from the standard scheme of 1/S- 
expansion H, [|, |£| since it gives a possibility to take into account exactly the kinematical 
relations between powers of spin operators on each site. We use the coherent state 
approach (see, e.g., Ref. 0) to write down the partition function in the form 



Z = J Direxp j iS J dr(l - cos^)|^ - (tt\H\tt) 



(9) 



where 7r is the unit-length vector, $ and <p are its polar and azimuthal angles respectively, 
1 7r) = exp(i$S y + itpS z )\S) are the coherent states (S Z \S) = S\S)). To construct the 
l/S'-expansion we rotate the coordinate system around the y-axis through the angle 9. 
The Hamiltonian ([!]) takes the form 



U 



— ^2 SjSj 

Z (ij) 



(10) 



Lm m'=—l 



vv v Rm (l + m)\(l-\m'\)\ A\ t 

I \ n _M/; , I /IM .\m'\ a mm 



,{9){0 



\m'\\ 

i H 



\ (I - m)\{l + \m'\)\ A 1 ™' 1 

where d l mm , (9) are the Wigner matrices of the rotation group irreducible representation, 

(l-m)\ 1 



AT 



([m] 



(/ + [m 
m for m even and [m 



ii 



m + 1 for m odd), for I < 4 we have K\ 



and the tilde sign here and hereafter is referred to the rotated coordinate system. 
Since the partition function (^) is invariant under rotation of the states \ir), it is 
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convenient to use the coherent states defined in the same coordinate system, i.e., 
1 7?) = exp(i$S y + i<pS z )\S) with S Z \S) = S\S). The advantage of using the coherent 
states is the simple form of the averages of the tensor operators (|2|) over \tt). By direct 
calculation we obtain 

(7r|0, ro |7f) = S l A?If l {co8'&) cosnvp (12) 

where Pj n {x) are the associated Legendre polynomials, the factors Si = S(S — 1/2)... [S — 
(I — l)/2] take into account properly the kinematical relations on each site. In 
particular, the second-order anisotropy term vanishes for S = 1/2, and the fourth- 
order for S = 1/2, 1, 3/2, as it should be (unlike the results of boson representations in 
Refs.p, ^ |6|). Using ([12]) we obtain for the case B™ = Bf5 m o under consideration the 
result 

JC2 

(mm= --|-E*i*i ( 13 ) 

+ EE E ^^^^pl m l(cos0)pH(cos^)cosmy, 

Further calculations are performed in the same line as in Ref. J7|. Representing 
cos^ = vT~— sin 2 § and expanding in sin^ we obtain the l/S'-expansion of the partition 
function. It should be stressed that we retain the factors Si, as well as S'-dependences 
in (f|), non-expanded. By performing decouplings, terms of third order are reduced to 
linear ones, and terms of fourth order to quadratic ones. The requirement of absence of 
sin ^-linear terms leads to the result for the cone angle 9 



2 n ^ 
cos 9 = - 
7 



l-X + Y--^(l- — + 6X + Y 
10 BlSi V 2S 



(14) 



where 



X = (t^) (sin 2 cos 2 <P) = T, %^ coth ^ 



Y = (vrj) = (sin 2 * sin 2 <p) = £ J ° ^oth ^, (15) 

and the "bare" magnon spectrum reads 



£ q = S J (J - J q )(J - J q + A /5) (16) 



Ar 



3B%S 2 cos 29 - 10B° 4 S 4 {28 cos 4 9 - 27 cos 2 9 + 3) 



A being the energy gap. The corrections in (jT3|) can be collected into powers in the 
same way as in Refs. fli~0| , [TT|j to obtain the correct description of thermodynamics at not 
too low temperatures. (In the presence of higher-order anisotropy this is essential since 



6 



the coefficients at X, Y increase as ~ Z 2 /2 with anisotropy order.) Then we have 



2 D ^ Zx 

cos e = -— 



Y 



1 Bl{T)S 2 
10 Bl{T)S A 



where 

B° 2 (T) = Z X Z Y B° 2 , Bl{T) = Z\Z Y B\ 
are the temperature-renormalized anisotropy constants, 

Z x = l + — - X, Z Y = 1 + — - Y 
2S ' 2S 



;i7) 



;is) 



(19) 



The relations (H) extend the results of Refs. (@, |, g |) to the case where spins make 
a non-zero angle with the z-axis. The renormalized gap has the form 

A = 6 cos 29 B°S 2 - 205°S 4 [3(1 - 7X) - 3 cos 2 9 (9 - 56X - 7F 
+28 cos 4 9 (l-QX-Y 



196 sin 2 cos 2 9 



xE 



3E 2 g 2 (Jo - Jk) ~ lQ^(g4/^)A cos 2 g 
ivl + S^Jo- J k )(J - Jk + Ao/5) 



(20) 



where X = X — 1/(25*), K = K — 1/(2S). After introducing the temperature- 
renormalized second- and fourth order anisotropy parameters D(T) and D'(T), 

D{T)S 2 = 3B° 2 {T)S 2 - 30E°(T)S 4 , 
D'(T)S 4 = 35B° i (T)S 4 (Z Y /Z x ), 



(21) 



the expression for cos# coincides with that of the phenomenological theory (P). 
Collecting again corrections in (|20"D into powers, we obtain for the renormalized gap 
in the notations ( pT|) the expression 



- 2 

!2) 



A = 2D{T)S 2 cos 29 + 2D'{T)S 4 cos 4 + 6D'(T)S 4 sin 2 cos 2 



196 sin 2 9 cos 2 J2 



3B^S 2 (J - J k ) - 10ff2(S 4 /S)Aocos s 
u; 2 + S 2 (J - J k )(Jo- Jk + Ao/S-) 



which also coincides with that obtained in the phenomenological theory except for the 
last term. Note that at 9 > the renormalizations (21) are present even at T = 0, 
which should be taken into account when treating experimental data. 



4. Orientational phase transitions 

Now we can pass to description of possible orientational phase transitions. Consider 
first the case of a small enough constant B® (or, equivalently, D'), so that cos 2 6^0) is 
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close to unity. Then cos 2 8{T) increases with temperature and there occurs a transition 
to the easy-axis phase at the point determined by 

3Zx 



7 Z 



Y 



i 1 B°(T)S 2 
10 Bl{T)S A 



(23) 



In the opposite case of a large enough B%, cos 2 6 l (0) is small, and cos 2 8{T) decreases 
with temperature, so that at the point where 

B° 2 (T)S 2 = 105° (T)^ (24) 

a phase transition to the easy-plane phase occurs. Thus one can expect that there exists 
the critical value 8 C : for 8q = 8(0) < 8 C we have a decrease of 8(T) with T and the phase 
transition from cone to easy-axis phase, while for tt/2 > 8 > 8 C we have an increase 
of 8(T) with T and the phase transition from cone to easy-plane phase. The numerical 
computations for the simple cubic lattice (see Fig.l) yield 8 C ~ 50°. Fig. 2 shows the 
corresponding temperature dependences of the anisotropy constants D(T), D'(T). For 
simplicity, J q is taken for the simple cubic lattice. 



The phase transitions described by Eqs. (^) and (p4|) are analogous to 
those in the phenomenological theory of Ref. |l[ that occur at D(T) = and 
D(T) = —2D'(T)S 2 , respectively. However, unlike the phenomenological approach, 
microscopical consideration leads to either cone to easy-axis or cone to easy-plane 
transition with increasing temperature, depending on the zero-temperature value of 
8. At the same time, the transition from the easy-plane to easy-axis structure (through 
the intermediate cone phase) cannot be explained by purely magnetic renormalizations 
of anisotropy constants. 

The result flT4] ) gives the mean- field values of the critical exponents (e.g., (3 = 1/2) 
for both the ground-state and temperature orientational phase transitions. Unlike the 
systems discussed in Ref. [0, the system under consideration has the dynamical critical 
exponent z = 2 (i.e., single excitation mode with nearly quadratic dispersion is present). 
Thus the upper critical dimensionality for the ground-state QPT is d^ = 4 — z = 2. In 



this respect, the system is analogous to XY model in the transverse magnetic field [12 
A characteristic feature of such systems is the mean-field behavior of critical exponents 
both above and below the critical dimensionality. For (hypothetical) systems with d = 2 
logarithmic corrections to ground-state properties near QPT are present (see, e.g., Ref. 
13fl). At the same time, the upper critical dimensionality for the temperature phase 



transition is d^ T = 4, and at d < d^ T the temperature-transition critical exponents differ 
from their mean-field values. 

Now we discuss the experimental situation. In Gd (see, e.g., Refs. 0, [14|| ) 
the orientational phase transition from cone-phase to easy-axis phase is observed at 
T c = 240K. The temperature dependence of the cone angle at T < T c (and also of 
magnetic anisotropy constants) is non-monotonous, unlike the results obtained in Sect 
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2. This complicated situation is connected with the absence of orbital momentum and 
smallness of anisotropy in gadolinium. 

In holmium the low-temperature phase is conical spiral one, the angle of the cone 
changing from « 80° to 90° in the temperature interval — 20K. The spiral angle makes 
up about 30°. Since the sixth-order anisotropy is important, we use the Hamiltonian 



^Ho = n + B l E(0 6 °), + Bt £(0 6 6 ) 4 (25) 

i i 

The hep lattice is not of a Bravais type. However, if we neglect the optical mode (which 
is possible at T T/v = 133 K) one can put (see, e.g., Ref. M) 



J q = 2J 



cosq x + 2 cos(q x /2) cos(\/3g J/ /2^ 



+ 2 J cos — 



exp(iq y / y/3) + 2 cos -|- exp(— iq y /2\/Z) 



(26) 



The parameters of the Hamiltonian were taken from Ref. ||15|| : J = 0.65K, J' = 
0.6 J, B% = 0.35K, Bl = 0, B% = -1.1 • 10- 5 K, B$ = 1.07 • 10- 4 K (note that our 
value of B® includes also renormalization due to dipolar anisotropy). For simplicity, we 
restrict ourselves to a collinear magnetic structure (this is justified by that the spiral 
angle in the rare earths is small, especially at low temperatures). The calculations 
with the Hamiltonian ( p5| ) are completely analogous to those in the previous Section. 
Calculated dependence of the cone angle is compared with the result of the mean-field 
approximation and experimental data in Fig.3. One can see that our results improve 
somewhat those of the mean-field theory where the temperature dependence of the 
anisotropy constants is given by (§). 

To conclude, we have formulated a consistent spin-wave approach to description 
of thermodynamic properties of anisotropic magnets at low temperatures. The 
renormalizations of the anisotropy constants and spin-wave spectrum for an arbitrary 
cone angle are calculated. This gives a possibility to describe the orientational phase 
transition between the cone and plane phases. 

We are grateful to J.Jensen for comments concerning the experimental situation in 
holmium. 

Figure captions 

Fig.l. The theoretical temperature dependences of the cone angle 9(T) for S = 7/2 
and different values of second-order anisotropy: D/J = 0.004; 0.005; 0.006 from upper 
to lower curve. The value of D'/ J is 3.7 • 10 -4 . 

Fig. 2. The temperature dependences of the anisotropy constants D(T), D'(T) 
corresponding to Fig.l. 

Fig.3. Calculated dependences of the cone angle in the mean-field approximation 
(short-dashed line) and renormalized spin-wave theory (RSWT, long-dashed line) as 
compared with experimental points for holmium (Refs. |2|, 
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